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The singularity of quantum correlations, such as quan- 
tum entanglement(QE) and quantum discord(QD), has 
been widely regarded as a valuable indicator for quan- 
tum phase transition(QPT) in low-dimensional quantum 
systems. In this paper, for an L x 2 spin ladder system 
with ring exchange, we find that the singularity of QD 



(or QE) could not indicate the critical points of the sys- 
tem. Instead, the QD shows a novel odd-even effect in 
some phases, which can be used to detect the phase 
boundary points. The size effect is related to the sym- 
metry breaking of the ladder. 
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1 Introduction The relation between quantum corre- 
lation and quantum phase transition (QPT) has been exten- 
sively investigated in many -body quantum physics. IU|2]|3] 
|4ll51|gll7l[8ll9l[TOl[T2l[T3ll29l It is first found by Osterloh et 
al. that the quantum entanglement (QE) — a kind of quan- 
tum correlation — shows a singular point in the vicinity of 
the QPT point of the transverse-field Ising model. |3] Since 
then, the ability of QE in detecting QPT has been inves- 
tigated in many quantum systems. Il4ll5ll6ll7l[8ll9l l29l How- 
ever, QE is not the only nature of quantum correlation, 
thus Olliver et al. [ITO proposes to use the quantum dis- 
cord (QD) to quantify all the quantumness of correlation 
present in a quantum state. |fT0l[TT1[T2l[T3l[T4l[T5l[T6l[T7l 

Just like QE, QD has also been used to study the QPTs 
in several one-dimensional quantum spin chains, such as 
the transverse Ising chain and an anti-ferromagnetic XXZ 
chain. lfl0l[T2l[T3l In most cases, these two measures of 
quantum correlation shows similar behavior in the phase 
transition. For example, a discontinuity of QE or QD can 
be used to detect a first-order phase transition point, while 
a discontinuity or a divergence in the first derivative of QE 
or QD would be related to a second-order phase transition 
point (or just called the critical point). lfT2l Though QD and 
QE have a similar capacity in detecting QPT, they have a 



fundamental difference, that is, QD can survive in separa- 
ble states.) 14] As a result, QD can capture the signal of 
QPT in separable states while QE cannot. For example, 
in an infinite XYZ model, the two-spin state of the sys- 
tem is separable thus pairwise QE vanishes, however, QD 
survives and indicates the QPT point of the system very 
well.lGJO) 

In addition to these simple models, recently novel 
quantum phases in complex spin systems have been stud- 
ied extensively, such as models with multiple-spin ex- 
change interaction, f 21} It's believed that multiple-spin 
exchange interaction plays an important role in under- 
standing the magnetic properties in several materials such 
as two-dimensional compound La9CttOzi. ll20l Among 
these models, the 5=1/2 spin ladder with four-spin ring 
exchange(the topology is shown in Fig. 1) has attracted 
special attention. I22ll231l24ll25ll26ll27ll Firstly, the two-leg 
ladder can be regarded as an intermediate topology be- 
tween one- and two-dimensional lattices, thus can give 
insight into the behavior of two-dimensional systems. Sec- 
ondly, due to the four-spin ring exchange interaction, the 
system has a very rich phase diagram. Using various cor- 
relation functions as order parameters, Lauchli et al. have 
found six phases, and the phase diagram of the system is 
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illustrated in Fig. 2. The non-symmetry-breaking regions 
contains four phases, i.e., a ferromagnetic (FM) phase, a 
rung singlet (RS) phase, a dominant vector chirality (VC) 
phase, and a dominant collinear spin correlation (CSC) 
phase, while the symmetry-breaking regions contains two 
phases, i.e., a staggered dimers (SD) phase and a scalar 
chirality (SC) phase. 12T1I251 We mark the two first-order 
QPT points as B\ and 62, the three critical points as #3, 84 
and #5, and the crossover point as 9q. 

Recently the QPTs of the ladder have been studied 
through quantum information theory. [22. 27. 28] The first- 
order QPT points 9\ and #2 are identified ambitiously by 
the sudden change of the entanglement concurrence. El 
The critical point 64, which is a highly symmetric point of 
the system, is identified by the size-independent extremal 
point of entanglement entropy. |22j For the crossover point 
9q, recently we have provided a very effective approach to 
identify its location by analyzing the first-excited state of 
the system. [28 1 Thus, in this paper, we will just pay our 
attention to the other two critical points #3 and 65. 

To investigate the entanglement properties of the lad- 
der, Song et al. have calculated the entanglement concur- 
rence |fT8l for several two-spin subsystems and the entan- 
glement entropy (9| for different block geometry, and the 
relation between the QE and the phase diagram is dis- 
cussed in depth. 1 22 1 However, neither the concurrence nor 
the entanglement entropy shows any singularity or extreme 
point at the phase boundary (#5) between the SC phase and 
the VC phase. Especially, in most areas of the SC phase and 
VC phase, the two-spin states of different subsystems (two 
spins on the rung, leg, and diagonal bond) in the ladder 
are all found to be separable states, thus the concurrence 
vanishes and shows no signal for 85 at all. The situation 
is somewhat similar to the above mentioned infinite XYZ 
model. Thus, it is natural to ask: would the QD capture 
the signal of the phase transition in these separable states'? 
Moreover, as far as we know, the existing works on QD are 
mainly limited to very simple models, and the connections 
between QD and novel phases in complex systems are still 
less well understood. Thus, it would be valuable to study 
the properties of the QD in the ladder. 

In this paper, firstly we give a brief introduction of QD 
in Sec. 2, then we describe the Hamiltonian of the ladder 
and some numerical details in Sec. 3. The main results are 
shown in Sec. 4 and some discussions are given in Sec. 5. 

2 Quantum discord QD aims at characterizing all 
the quantumness in a quantum state. 11111141 Its definition 
is based on the two quantum versions of the classical mu- 
tual information. (T9) For a classical system (or a state) 
AB, the total correlation between the subsystems A and 
B can be expressed as Ia.b = Ha + Hb — Bab or 
Ja,b = Ha — Ha\b, where Ba (Bb, Bab) is the Shan- 
non entropy, and Ba\b is the conditional entropy. One can 
prove that Ia.b and Ja.b are equal to each other. Now let's 
extend the definition of Ia.b and Ja,b to a quantum state 
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Figure 1 (Color online) The topological structure of the 
two-leg ladder. 



described by the density matrix pab- For Ia.b, by replac- 
ing the Shannon entropy with the von Neumann entropy 
S(p), one can easily obtain the quantum mutual informa- 
tion x(pab) asunnnnnsna 



Z(pab) = S(p A ) + S(p B ) - S(pab)- 



(1) 



The quantum generalization of Ja,b is slightly diffi- 
cult because of the conditional entropy Ba\b, thus, firstly 
the quantum conditional entropy S(p\{Bk}) — the quan- 
tum version of Ba\b — has to be defined, where {B^} 
is just a complete set of projectors. 111511161 Then a vari- 
ant of quantum mutual information can be defined as 
JjB k} {pAB) = S(p A ) - S(p\{B k }). Following Ref. 
[11|, the quantity 



J(PAB) ■= Sup J{B k }(PAB) 
{B k } 



(2) 



is interpreted as a measure of classical correlation. Finally, 
QD is just defined as the difference between the quantum 
mutual information X(pab) and the classical correlation 
J(Pab)- 



Q(pab) = X(pab) - J(pab)- 



(3) 



For a state containing quantum correlation, QD is gen- 
erally non-zero. However, if pab reduces into a classical 
state, one can prove that I(pab) and J{pab) would just 
reduce to Ia.b and Ja.b, respectively, thus the QD van- 
ishes. 

For a general state, the above mathematical definition 
of QD is somewhat complex, with a numerical optimiza- 
tion procedure involved. However, for two-qubit states ex- 
pressed in the following form, 
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the QD can be simply expressed as spin-spin correlation 
functions, and the numerical optimization procedure can 
be omitted. [ 16| Because of the symmetry of the ladder con- 
sidered in this paper, the two-spin reduced density matrix 
takes the form in Eq. |4] thus our calculation follows the 
formula in Ref. Ifl6l . 

3 Spin ladder with ring exchange We consider the 
following SU(2) invariant Hamiltonian defined on a L x 2 
two-leg ladder as lT2D 



+^Eti(4,+i + 



i+l -r 
P- 1 

i, i + l 



Sd,i • Sd,i+i) (5) 



where the two legs of the ladder are marked as u and d, and 
the rungs are labeled as i = 1, • • • , L(see Fig. 1). The first 
term in the above summation describes the interaction on 
the rungs, the second term is the interaction along the legs, 
the third term is the four-spin exchange interaction in the 
plaquettes of the ladder, and J r , Ji and K are the corre- 
sponding coupling constants. S u ^{Sd.i) is the spin-1/2 op- 
erator defined on the z-th spin in leg u(d). Pi i i+i(P i ~^ 1 ) is 
the four-spin ring operator, which rotates the four spins in 
the i-th plaquette clockwise (counterclockwise), i.e., J 
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(6) 



In this paper the periodic boundary condition is imposed by 
setting S u i = S U: l and Sd,i = Sd,L- In addition, follow- 
ing Ref. l2TI . we set J r = Ji = cos(#) and K = sin(9). 
The phase diagram is shown in Fig. 2. 

We will investigate the behavior of quantum correla- 
tion in the vicinity of 6*3 and #5. For finite-size ladders with 
L = 6, 8, 10, 12, firstly, we use the exact diagonalization 
method to obtain the ground-state wavefunction. Then we 
calculate the reduced density matrix of the concerned sub- 
system by tracing out all other degrees of freedom. Finally, 
QD can be obtained analytically according to Ref. |[T6l . 

We have calculated the QD for the two spins on the 
rung(sites S u ,i and S d ,i), leg(sites S u ,i and S u ,i+i), and 
diagonal bond(sites S U) i and Sd,i+i) in the ladder, and the 
most interesting results are the behavior of QD on the rung. 
Thus, firstly, we limit ourselves to the two spins on the 
rung, and discussions about other two-spin subsystems will 
be given in Sec. 5. 

4 Odd-even effect between phase boundaries 9 3 

and 65 The QD between the two spins on the i-th rung of 
L x 2 ladders with L = 6, 8, 10, 12 is shown in Fig. 2. 

One can see that at the QPT points 6\, 82 and 84, the be- 
havior of QD is similar to the entanglement reported in pre- 
vious studies. Il22ll that is, QD is discontinuous at the first- 



oder QPT points 8\ and 82, and shows a size-independent 
extremal point at the highly symmetric point 64. 

Now we try to find the signal for #3 and 85. From Fig. 
2, it is expected that QD would not show any singularity 
or extremal point at these two points when L is increased, 
thus we need to change our strategy. We find that the finite- 
size effect in the symmetry-breaking regions (8 S [83,85]) 
is obviously larger than that in the non-symmetry -breaking 
regions (8 £ [83 ,85}), which may be useful to detect #3 and 
5 - 

In the non-symmetry-breaking regions, the size effect 
is very small. For the FM phase, for example, in the rig- 
orous bounds — 7r < 8 < — 5 , it is clearly that the fer- 
romagnetic state minimizes the energy on each plaquette 
locally, thus the finite-size effect is relatively small. For 
the RS phase, let's choose the point 8 = 0, i.e., the spin 
ladder with only anti-ferromagnetic bilinear couplings. It's 
well-known that the ground state can be well approximated 
by the product of local rung singlets, thus the size effect 
should be small in this phase. In the main region of the VC 
phase and the CSC phase, QD also shows very little size 
dependence. 

In the symmetry-breaking regions, peculiar size effect 
is observed, that is, the behavior of QD with L = M is 
notably different from that with L = 4t + 2. The physical 
origin of the peculiar size effect can be understood as fol- 
lows. In the SD and SC phases, the translational symmetry 
is broken and the "unit" of the ground state consists of two 
rungs, thus, the number of "unit" is even for L = 4i and 
odd for L = 4i + 2. Then, the wavefunction and quantities 
such as QD would show an odd-even dependence upon the 
number of "unit" of the ground state as follows: 



l* 9 > = { 



\ip L ) fori = Ai 
\tp L ) foiL = 4i + 2 



and 



Q 9 = { 



QtP l for L = 4i 
Q VL forL = 4i 



(7) 



(8) 



In order to see the size effect more clearly, for different 
phases we numerically calculate the QD for the two lowest- 
lying energy eigenstates, and several typical results are 
shown in Fig. 3. For a given 8, it is obvious that these (L,Q) 
points belong to two curves, thus we label one of the two 
curves with squares and dot line (i/jl) and the other with 
circles and dot line (<pl)- In addition, for the ground state, 
we mark the corresponding points (L,Q g ) by red triangles. 
In non-symmetry -breaking phases(see Fig. 3 (a)(e)(f)), one 
finds that the ground state is always \(Pl), while \4>l) is 
always an excited state. In symmetry-breaking phases(see 
Fig. 3 (b)(c)(d)), however, with the increase of L, (L,Q g ) 
jumps between the two curves, which means these two 
states take turns to serve as the ground state as the change 
of L, just as described by Eqs. |7]and|8] 
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Figure 2 (Color online) Ground-state QD of two spins on the i-th rung as a function of 8 in L x 2 spin ladders. The 
phase boundaries are marked as 8i. Arrows with solid line, dash line and dot line denote first-order QPT points {8\ and 
62), second-order QPT points (6*3, 84 and #5), and the crossover point (f9 6 ), respectively. The values of $i from Ref. ED 
are shown in the figure. 



This jumping behavior in Fig. 3 can be used to detect 
the phase boundary points f9 3 and 65 . However, much man- 
ual work would be required to obtain Fig. 3 for all the 9. 
Based on Fig. 2, there is an alternative approach to locate 
83 and #5 roughly but very efficiently. First, one defines 
a parameter q L (9) := \Ql{9) - Ql-2(8)\- From Fig. 2 
one sees that with the increase of L, qh{8) vanishes very 
rapidly in the RS and VC phases. For a given 9, with L 
large enough, suppose qi,{ff) is smaller than some thresh- 
old value 8, we say that the system is in the non-symmetry- 
breaking region, thus #3 and #5 can be identified approxi- 
mately. In fact, when using L = 12 and <5 ~ 10~ 4 , #3 and 
05 are found to be 0.077T and 0.397T, respectively, which is 
consistent with previous studies. lt2~D 

From Fig. 3 one can see that the odd-even behavior de- 
creases with increasing system size. Therefore, one might 
wonder, if this criterion is useful for large systems. Keep in 
mind is that, for a given 9, if the odd-even behavior is ob- 
served for a small L, the system should be in the symmetry- 
broken phase, even if the behavior becomes blurred when 
the size of the system is large. Thus, for a finite L, the odd- 
even behavior would be of practical use to detect, at least, 
the "inner boundary" of the symmetry-broken phases. In 
other words, the above values can always be regarded as 
an upper (lower) bound for the critical values #3 (#5) of the 
ladder systems. 

5 Discussions and summaries In many systems, 
the singularity of ground-state wavefunction is widely used 
to detect the QPTs (the singularity of quantum entangle- 
ment and quantum discord is also inherited from the wave- 
function). However, in this paper, in a finite-size ladder 
with ring exchange, the wavefunction keeps analytical at 
6*3 and 85. Fortunately, we find that the odd-even effect in- 
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Figure 3 (Color online) The QD of the two lowest-lying 
energy eigenstates of the system for several 8. The red tri- 
angle denotes the QD of the ground state. In the SD and SC 
phases, the ground-state QD jumps between two curves, 
which indicates that two states \4>l) and \<pl) take turns to 
serve as the ground state. 



duced by the symmetry breaking of the wavefunction can 
be used to detect the QPT point of the model. 

The method may be a little rough compared with other 
large-scale numerical methods dealing with the infinite- 
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size systems. II27II However, it captures the approximate lo- 
cations of the QPTs with the negligible amount of calcula- 
tion (the results are obtained using a standard PC in about 
two days). Note that a valuable estimation of the QPT point 
would always save much time before large-scale numerical 
calculations are adopted to figure out the accurate location 
of the QPTs. In addition, as we have shown, the size depen- 
dence of QD provides a special perspective to understand 
the physical nature of the symmetry-breaking phases. 

As QD and the concurrence are closely related quan- 
tities, we'd like to compare the behavior of QD and the 
concurrence in the system. Both of them describe the quan- 
tum correlation between the two spins in the system. QD 
is meaningful in the whole parameter space: it can be used 
for both entangled states and separable states. The concur- 
rence, as has been mentioned, would vanish in a separable 
two-spin state even if the whole system is in an entangled 
state. [ 10] That's why the QD can be used to detect #5 while 
the concurrence fails. EZ) Thus, when used as phase transi- 
tion indicator, QD has a broader scope of application than 
the concurrence. 

Our studies are limited to the QD of the rung. In fact, 
for other two-site subsystems, such as two spins on one leg 
or the diagonal bond in the ladder, the QD has also been 
calculated(the result has not been shown in this paper), 
however, QD does not show such obvious odd-even behav- 
iors, thus there is no clear signal for the boundaries #3 and 
65. The reason is as follows. The observed finite-size effect 
always contains two sources: source (i) is just the selection 
between \iJ)l) and |(/?i), described by Eq. [7] and source 
(ii) comes from the convergence speed of the wavefunc- 
tion I^lHI^l)) when increasing L. Source (i) just emerges 
in symmetry-breaking systems, while source (ii) is gener- 
ally present in low-dimensional quantum systems. To get 
a clear observation of odd-even effect from source (i), the 
effect from source (ii) need to be weak enough, that is, ob- 
servables (such as QD) calculated from |</?l)(|V'l)) need 
to converge very fast. Because of the special topology of 
the ladder, it is expected that the speed of convergence 
for inter-rung correlations would be slower than that of 
the intra-rung correlation. Il30l As a result, for QD of the 
two spins on the rung, the effect from source (ii) is weak 
enough thus the effect from source (i) can be revealed ap- 
parently, while for QD between different rungs, the effect 
from source (ii) is very large thus the effect from source 
(i) becomes blurred. Thus, it is expected that QD of the 
spins on the rung, rather than other two-spin subsystems, 
is most likely to reveal the odd-even effect from source (i) 
and show signals for #3 and 65. 

The results obtained in this paper stimulate us to study 
the size effect of quantum correlation in some other models 
with spatially inhomogeneous phases. OTll 
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